gsem (Tability -> a6@-1 a7 a8 a9 a10 , family(ordinal) link(probit) ) /// 
(a6 <-  agem male ,) ///
(a7 <-  agem male ,) ///
(a8 <-  agem male ,) ///
(a9 <-  agem male ,) ///
(a10 <-  agem male ,) , ///
gin(cons) latent(Tability)  startvalues(constantonly)  technique(nr) difficult

*predict tab if e(sample), latent(Tability)
*label var tab "empirical Bayes means for Tability: monthly task a6-a17"

gsem (Intability -> a18@-1 a11 a12 a15 a16 a17 a19 a20 a21 a22, family(ordinal) link(probit) ) /// 
(a11 <-  agem male ,) ///
(a12 <-  agem male ,) ///
(a15 <-  agem male ,) ///
(a16 <-  agem male ,) ///
(a17 <-  agem male ,) ///
(a18 <-  agem male ,) ///
(a19 <-  agem male ,) ///
(a20 <-  agem male ,) ///
(a21 <-  agem male ,) ///
(a22 <-  agem male ,) , ///
gin(cons) latent(Intability)  startvalues(constantonly)  technique(nr) difficult

*predict intab if e(sample), latent(Intability) intp(11)
*label var intab "empirical Bayes means for Intability: monthly task a18-a22"


gsem (Int_tc -> a23@-1 a13 a14 a24 a25 a26 a27, family(ordinal) link(probit) ) /// 
(a13 <-  agem male ,) ///
(a14 <-  agem male ,) ///
(a23 <-  agem male ,) ///
(a24 <-  agem male ,) ///
(a25 <-  agem male ,) ///
(a26 <-  agem male ,) ///
(a27 <-  agem male ,) , ///
gin(cons) latent(Int_tc)  startvalues(constantonly)  technique(nr) difficult

*predict int_tc if e(sample), latent(Int_tc) intp(7) tol(1e-6)
*label var int_tc "empirical Bayes means for int_teacher child: monthly task a23-a27"
